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ABSTRACT: First-principles calculations combined with XRD simulations are performed to P™r0P 
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systematically study crystal structures, bonding characteristics and electronic structures of **@. i s rem 
Li,CoSiO, (x = 2.0, 1.5, 1.0) polymorphs with symmetries Pmn2,-DP and Pbn2,. The ‘cot è ae | 
calculated average voltages by lithium extraction agree well with available experiments. CoO, : = Ea B y l 
tetrahedron is the key structural unit to track the process of delithiation. The oxidation of w A 5 Ey 
CoO, tetrahedron results in a special pattern of bonding characteristic, which corresponds to i bee 3% Serre 
spin ordering and may be observable in XRD spectra according to simulation. We find Pom O E] | z 
delithiated phases are intrinsic Mott insulators, electronic band gaps change from Mott— o2@ Q 
Hubbard-type to charge-transfer-type during lithium removing. The swapping of near-gap wet 4 lk 
states is associated with the contraction of the oxidized CoO, units, indicating Peierls a . | 
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distortions that may be the physical origin of capacity degrading of Co—silicate chemistry. i i T arid 


1. INTRODUCTION 


The urgent need to double energy density from the current 
state-of-the-art lithium-ion batteries has drawn intense interest 
in finding or designing a new chemistry beyond oxide cathode 
materials such as LiCoO, and LiMn,O,.' ° High capacity in 
the desired voltage range of electrolytes with better safety and 
longer cycle life is continuously the preliminary requirement of 
the new chemistry for application to scalable energy storage in 
electric transportation and power-grid.*~° Polyoxyanion 
compounds, typically, LifePO, proposed by Goodenough, et 
al, have received wide recognition as safe cathode materials.’ 
Extensive studies on olivine phosphate chemistry have been 
carried out in order to understand the physical origin of 
electrochemical performance of polyoxyanion frameworks.*~"* 

Transition metal silicates Li,MSiO, (M = Fe, Mn, Co, Ni) 
have attracted significant attention as a new type of 
polyoxyanion cathodes since 2005.'*~'? All cations in 
Li,MSiO, form tetrahedral units (XO,)” (X = Li, M, Si), 
very different from octahedral units in olivine phosphate 
chemistry. All silicate formulas indicate high theoretical 
capacities over 300 mAh/g if two lithium ions extracted. Fe— 
and Mn-silicates have achieved nearly 2 Li* extraction in 
recent works, 7°? thanks to continuing advances in synthesis of 
nanostructure and carbon-coating.””” However, electrochem- 
ical performance of Fe— and Mn-silicates seems unstable and 
their working voltages are low, under 3 V, not preferred with 
current electrolytes for high energy density. Recent calculations 
have carefully analyzed the bonding characteristic of Li,FeSiO, 
in order to determine the critical structural units that impact 
electrochemistry of tetrahedral framework.'* It has been 
proposed to modify silicate chemistry through synergy effects 
as did the transition metal ternary oxides.'® Such an attempt in 
silicates was mostly limited to the mixing of Fe and Mn, and no 
improvement was reported in literature so far.” 


It is worth pointing out that Co also play a key role in the 
solid solution ternary oxides, Li(NiCoMn),/30, and Li- 
(Nig. gCo9,1sAlg9s)O2. We notice that Co silicate could 
discharge at around 4 v,?3 almost 1 V higher than both 
Fe— and Mn-silicates and as-prepared Co—silicate polymorphs 
did not present voltage degradation as always did Fe— and 
Mn-silicates.”*”> However, Co-silicate chemistry is still 
struggling to deliver high capacity. Therefore, it is important 
to understand the intrinsic electrochemistry of Co—silicate and 
its difference from other silicate chemistry. Previous first- 
principles calculation indicated a high average voltage at 5.0 V 
would be required to extract fully two lithium ions and average 
voltages for partially delithiated phases distributed in a wide 
range for Co-silicate.°°-** Those predictions disagreed with 
recent experiments on carefully characterized Co—silicates with 
symmetries Py-Pmn2,, PrPbn2, and Yor P2,/n P25 In our 
recent work, we identified all the three synthesized polymorphs 
were in 3D framework lattices, containing no 2D layers existing 
in Fe— and Mn-silicates, whose crystal structures have been 
used as the starting template in previous Co-—silicate 
calculations. And the difference of lattice dimensions in Co— 
silicate polymorphs can be well distinguished by their X-ray 
diffraction (XRD) spectra.” 

This work, mainly by using first-principles calculations and 
XRD simulations, aims to offer further details on the effect of 
the crystal structure and electronic structure of Co—silicate 
Li,CoSiO,(x« = 2, 1.5, 1.0) polymorphs on their electrochemical 
properties. The average voltages are calculated to compare with 
the available experiment. We find that CoO, tetrahedral arrays 
are a critical structure unit in the process of lithium 
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Table 1. Equilibrium Volumes, Total Energies, Magnetic Moments per Co, and Average Voltages for Li,CoSiO, and Delithiated 


Li,CoSiO, (x = 1.5, 1.0)* 


Li,CoSiO, 
Pmn2,-DP Pbn2, 
Q(A?)/formula unit 85.095 85.071 
Energy(eV) /formula unit —51.545 —51.547 
m (up)/Co ion 2.759 2.761 
Voltage-Cal.(V) 
Exp. (V)” 


Others-Cal.(V) 


Li, <CoSiO, LiCoSiO, 
Pmn2,-DP Pbn2, Pmn2,-DP Pbn2, 
86.452 85.093 87.302 87.378 
—48.564 —48.650 —45.559 —45.539 
2.747, 3.103 2.748, 3.112 3.166 3.166 
4.062 3.896 4.110 4.320 
3.7—4.0 3.5—4.3 4.0—4.3 4.3—4.6 


4.4%, 3.44, 4.135 4 


“Average voltages for the first 0.5 Li (per formula unit) extraction are from Li,CoSiO, to Li, ;CoSiO, and the second 0.5 Li (per formula unit) from 
Li, ;CoSiO, to LiCoSiO,, respectively. All calculations are in FM configurations by full GGA+U calculation including structural relaxing. ’Estimated 
from ref 23. “Reference 26. “Reference 24. “Reference 28. Reference 42. All other calculations were based on layered Pmn2, structures for one 


lithium extraction from Li,CoSiO, to LiCoSiO,. 


intercalation/extraction. Both atom magnetic moments and 
local lattice structures of CoO, units are investigated for 
delithiated Li,CoSiO,(« = 1.5, 1.0). The electronic band 
structures of delithiated phases are found to be intrinsic Mott 
insulators, changing from Mott—Hubbard-type to charge- 
transfer-type during delithiation. This may provide physical 
insights on the limited capacity of Co—silicates. We believe that 
it is essential to identify the key structural unit that is important 
to polymorph structure and redox potential for silicate 
chemistry in general. The understanding of Co—silicate may 
benefit not only its own electrochemistry but also its use as a 
potential component in solid solution silicates. 


2. COMPUTATIONAL DETAILS 


We have determined the key crystal structures of polymorphs: 
By-Pmn2,-DP, f,-Pbn2,, and 7o-P2,/n in our recent work via a 
detailed comparison with the carefully characterized exper- 
imental structures.” We have attached DP to Py-Pmn2, to 
explicitly indicate its XRD feature and its difference from 
layered Pmn2, structures. The primary cell contains 2 formula 
units for Pmn2,-DP symmetry and 4 formula units for both 
Pbn2, and P2,/n symmetries. All the three polymorphs are 3D 
framework lattices, with no 1D lithium ion channel as in 
LiFePO, or 2D lithium ion layers as in LiCoO, or Fe— and 
Mn-silicates. All of them can be distinguished by characters in 
XRD, Pmn2,-DP structure shows a double peak around 25° in 
20 while both £,-Pbn2,, and 7-P2,/n show similar one-triple 
peak patterns in similar angle ranges. Detailed structure 
information and comparison with experiment spectra for the 
fully lithiated Li,CoSiO, has been reported in our recent 
work.” Therefore, this work will focus on partially delithiated 
Pmn2,-DP and /,-Pbn2, to address the key relationship 
between polymorphic structure and electrochemical properties. 

All the calculations were performed using the Vienna ab 
initio simulation package (VASP)?! within the framework of 
density functional theory (DFT), employing the projector- 
augmented wave (PAW) method with the spin-polarized 
generalized gradient approximation (GGA) of Perdew— 
Burke—Ernzerhof (PBE), $’ °4 and adding the Hubbard 
parameter correction (GGA+U) to address the strong onsite 
Coulomb interaction for transition metals. DFT+U method and 
hybrid functional method are two popular schemes treating the 
self-interaction error in the conventional DFT framework. But 
DFT+U method is superior to hybrid functional method in 
predicting Li intercalation potentials, magnetic moments, and 
charge localization for all the positive electrode materials 
containing 3d electrons in transition metals, due to hybrid 


functional method being sensitive to the bonding environment 
of the oxygen, in a way not found in GGA+U.*” Extensive 
calculations have showed that GGA+U corrections can give 
quantitatively good agreement with experimental average 
voltages only through careful adjustment of values of U, 
which strongly depend on the elements and the valence states 
in materials.” For GGA+U calculations, the “effective” value 
(U-J) was set to be 5.0 eV for Li,CoSiO, (x = 2.0, 1.0), such 
parameters setting has been considered suitable in previous 
first-principles studies***°*” While for Li, ;CoSiO,, we found 
an effective value of 4.5 eV better for describing the formal 
charge separation between Co”* and Co**, due to the slight 
difference of localized d-electrons associated with Co** and 
Co** ions. A 500 eV cut off for the plane-wave basis set was 
used. Integrations in the first Brillouin zone were performed 
using Monkhorst-Pack special k-point mesh of 10 X 12 X 12 
and 10 X 6 X 12 for the primary unit cells of Li,CoSiO, 
polymorphs with symmetries Pmn2, and Pbn2,, respectively. 
Structural relaxations were performed with the total energy 
converged to 107° eV or better and all forces acting on ions 
found to be 0.005 eV/A. The same FM configurations were 
used to calculate delithiated Li,CoSiO, as in Li,CoSiO,. 

We used the standard first-principles voltage calculation,**® 
which is derived from the calculated total energies of host 
compounds 
7 — _ BlbisCoSiO,) - B(LI,CoSiO,) = (« = y)u(Li) 


(x — y) 

(1) 
where E(Li,CoSiO,) is the total energy of the initial phase 
Li,CoSiO,, E(Li,CoSiO,) is the total energy of the following 
delithiated phase Li,CoSiO,, H(Li) is the total energy of BCC 
lithium metal, (x — y) is the fraction of extracted lithium ions 
per formula unit, and V is the theoretical average voltage for 
lithium concentration between x and y, generally used as the 
intrinsic voltage of electrode materials of Li-ion batteries and 
compared to experimental voltage measurements in literature. 
Abundant first-principles calculations over past decade have 
showed the accuracy and predictive power of this voltage 
expression (1), often within 0.1—0.2 V of experimental voltage 
profiles for many cathode materials.*”* For Li; sCoSiO,, we 
compared two nonequivalent lithium-vacancy configurations 
for Pmn2,-DP, and nine nonequivalent configurations for 
Pbn2,. For LiCoSiO,, we considered three nonequivalent 
distributions for Pmn2,-DP and 35 representative patterns for 
Pbn2,. The lithium-vacancy configuration having the lowest 
energy is used to calculate the voltage profile by eq 1, and then 


corresponding delithiated structure is further analyzed to offer 
insights on the physical origin of electrochemical properties of 
Co silicates. 


3. RESULTS AND DISCUSSION 


3.1. Voltage Profile of Li,CoSiO, Polymorphs. Table 1 
shows that the voltage of Pmn2,-DP gently increases from 
4.062 to 4.110 V when extracting one full lithium ion per 
formula unit. On the other hand, the Pbn2, structure, while 
starting from a slightly lower voltage at 3.896 V for the first half 
Li’ extraction, jumps to 4.320 V during the second half Li* 
extracting from Li,,;CoSiO, to LiCoSiO,. This notable 
difference of voltage platform changes agrees with the 
experimental observation.” Up to cutoff 4.3 V, about one 
lithium per formula unit was extracted from Pmn2,-DP but only 
about 0.5 lithium extraction per formula unit for Pbn2,.”° To 
remove one lithium ion from Pbn2,, cutoff voltage had to be 
increased to 4.6 V. This confirms a higher voltage platform 
needed to extract the second half lithium ion from Pbn2,. 

Table 1 also shows previous voltage calculations of partially 
delithiated phases are distributed in a wide range for Co— 
silicates, 482342 ranging from 3.4 to 4.4 V and disagreeing with 
experimental trends. We notice that previous Co—silicate 
calculations were often using Fe— and Mn-silicates as the 
starting template, which contains 2D lithium layers.**°-*** 
As we have pointed out that there exists no such lithium layers 
in synthesized Co-silicate polymorphs according to XRD 
analysis.°° Voltage calculations indicate that the electrochemical 
properties of silicates also strongly depend on detailed 
tetrahedral networks.** The dimension of tetrahedral networks 
is the basic structural difference between Co-—silicates and other 
silicates in synthesized polymorphs. We have found that CoO, 
arrays are a suitable indicator of cationic ordering to track 
polymorphic structures of Co-silicates in XRD.’ Our 
calculations of voltage profiles indicate cationic orderings also 
remarkably affect details of voltage platforms. And our 
calculation suggests that Pmn2,-DP is within the fair voltage 
window of current electrolytes. 

3.2. Spin Ordering Formation with Oxidation of CoO, 
Arrays. Table 1 gives magnetic moments of 2.759 and 2.761 jug 
per Co in Li,CoSiO, for Pmn2,-DP and Pbn2,, respectively. A 
clear pattern shows up in magnetic moments with the change of 
valence states of Co ions during delithiation. All polymorphs 
confirm the formal valence Co”* with electronic configuration 
of 3d’ at a high-spin state, close to 2.78 up per Co in the same 
valence Co™* at the high-spin state of the olivine structure 
LiCoPO,.” However, the Co’ ion in Li,CoSiO, is in 
tetrahedral coordination with surrounding weak oxygen ligand 
field. The quasi-tetrahedral crystal field in Li,CoSiO, splits the 
5-fold degenerated states of the Co ion into two subsets, two 
low-lying e states and three high-lying t, states, different from 
the order (tag, ep) in the strong octahedral ligand field of the 
olivine structures. 

The redox pair Co**/Co** is conceptually believed important 
to electrochemistry of Co—silicates. As pointed out in a recent 
work,** magnetic orders are measurably preferred to the 
concept of charge-on-ions for the formal valences in the solid 
state. Table 1 confirms such a clear correlation between the 
spin states and the formal valences of Co ions in the course of 
oxidation. It is of interest that the magnetic moment 3.166 pg 
of the tetrahedral Co** in LiCoSiO, is again comparable to 3.27 
Hp of the octahedral Co** in olivine phosphate.” Delithiated 
Li, ;CoSiO, consists of half Co in the formal double and triple 


valences, respectively, as manifested by two types of magnetic 
moments, or formally called the charge separation of Co ions 
on the CoO, arrays. Parts a and b of Figure 1 show a detail 
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Figure 1. Layout of the spin ordering in delithiated Li, ;CoSiO,: (a) 
Pmn2,-DP; (b) Pbn2,. Green, blue, gray, and red balls correspond to 
Li, Co, Si, and O atoms with big balls on top layer and small balls on 
low layer for comparison, with not all ions shown. Magnetic moments 
are showed in units of 4p. 


pattern of the spin ordering (or magnetic ordering) for Sy and 
Pı Li, ;CoSiO,, respectively. In Figure 1a, different spin Co ions 
settle on different layers; in Figure 1b, both lay on the same 
layers. We find that the formal valence agrees with the local 
oxidation situation of CoO, arrays. The Co** ions with 3.1 4g 
are nearby the Li-vacancy while the Co”* ions with 2.7 jug locate 
distant from the Li-vacancy, confirming the formal charge of Co 
ions. Correspondingly, the magnetic moments of O ions split 
according to the oxidation states of nearby Co ions: in the 
range of 0.03—0.07 jig if near Co” ions and increasing to 0.10— 
0.20 lp if near Co** ions. Thus, tetrahedral CoO, as a group is 
better to track the oxidation/reduction process than does single 
Co ion. 

Table 1 shows that the volume difference between Jy and fy 
become large at Li,;CoSiO, but become small at LiCoSiO,. 
Pmn2,-DP presents a smooth volume expansion over the entire 
extraction while Pbn2, shows a steep expansion from 
Li, ;CoSiO, to LiCoSiO,, which also accompanies the voltage 
jump. This is correlating to the different pattern of charge 
separation that increases Coulomb repulsions among ionic 
tetrahedrons by the oxidation of Co% to Co** ions. Table 2 
reveals a clear evolution of bonding characteristic over the 
entire oxidation from Li,CoSiO, to LiCoSiO,. Overall, the 
average Si-O bond lengths and SiO, tetrahedral volumes 
remain almost unchanged, indicating a strong covalent bonding 
between Si and O atoms that plays a central role to keep the 
structural skeleton stable during lithium extraction. The average 
Li—O bond lengths show a slight change of about 1% and LiO, 
tetrahedral volumes change less than 2.5%. 

Table 2 shows a 20% volume difference between Co**O, and 
Co**O,, in according with about 7% difference in average Co— 
O bond lengths. The oxidation of Co**O, to Co**O, units, 
however, decreases electronic repulsions between Co and the 
surrounding O ions, mostly due to the transfer of electrons on 
one e orbitals; consequently, the bond length between Co and 
O sites becomes reduced. The Co ions in the CoO, units with 
shorter average Co—O bond lengths have higher magnetic 


Table 2. Variations of Average Bond Lengths (Li—O, Co—O, 
Si—O) and Tetrahedral Volumes (LiO,, CoO, SiO,) upon 
One Lithium Extraction from Li,CoSiO,, through 

Li, ;CoSiO,, to LiCoSiO, for Pmn2,-DP and Pbn2, 
Structures 


Pmn2,-DP Pbn2, 
Li,CoSiO, x=2 x=l5 x=1 x=2 x=l15 x=1 
dco-o (A) 1.995 1.993 1.870 1.997 2.006 1.871 
1.870 1.872 
dsi-o (A) 1.657 1.655 1.648 1.657 1.654 1.649 
di-o (A) 2.004 2.026 2.023 2.003 2.005 2.023 
Vcoo, (A®) 4.064 4.022 3.296 4.075 4.127 3.328 
3.294 3.302 
Vsio, (A®) 2.332 2.320 2.295 2.332 2.321 2.296 
Viio, (A®) 4.094 4.190 4.116 4.086 4.064 4.034 


moments at about 3.1 jg, while those in the CoO, units with 
longer average Co—O bond lengths have lower magnetic 
moments at about 2.7 jg (cf. Figure 1), in compliance with the 
inverse dependence of valence states of transitional metal M on 
average M-O bond lengths.*° In Pmn2,-DP, the Co**O, unit in 
Li, ;CoSiO, contract slightly from that of Li,CoSiO, while it is 
a slight expansion in Pbn2,. 

There is a concern on whether or not such a theoretical 
picture of spin ordering (or charge separation) of CoO, arrays 
can be observed.** Figure 2 predicts the bonding characteristic 
associated with the oxidation indeed produces visible features in 
simulated XRD spectra. An obvious splitting of the main peaks 
shows up in the spectra of Li,;CoSiO,, and due to the 
significant contraction of tetrahedral Co**O, units, a set of 
peaks between 25° and 30° in Figure 2 disappear or combine 
into a new one around 27° at LiCoSiO,. Such a spectral 
observance should confirm the fine effect of oxidation on the 
CoO, tetrahedral bonding network. 

3.3. Effects of Electronic Localization. Expanding lattices 
with shrinking CoO, units would prompt one to expect a large 
space for lithium ion movement, which might suggest a better 
rate capability and easy discharging for lithium move-in. 
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However, all synthesized polymorphs of Li,CoSiO, presented 
awkward discharging capacity even right after the first 
charging.” Hard electron transport has been blamed for the 
bad discharging. However, the nature of the electronic 
structures at delithiated phases (x = 1.5, 1.0) and the influence 
of the lattice change on the electron configuration upon the 
CoO, oxidation have not been addressed in details so far. 
Our GGA+U calculations find both delithiated phases (« = 
1.5, 1.0) are Mott insulators for both Jy- and £,-polymorphs 
and the nature of insulating gaps strongly depends on the 
oxidation states of Co ions. Parts a and b of Figure 3 show that 
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Figure 3. The partial density of states of Co ions on different 
tetrahedral CoO, sites (c. f. Figure 1) in delithiated Li, ;CoSiO,: (a) 
Pmn2,-DP; (b) Pbn2,. Fermi levels are at zero. 


gap states are similar for Jy- and /,—Li,;CoSiO,4. Pbn2, has a 
slightly bigger gap than Pmn2,. The density of states (DOS) of 
Co% in the upper panel of Figure 3 indicates a Mott—Hubbard- 
type insulating gap between sharp e- and t,-bands for the 
electronic configuration e“t,>. The occupied spin-up bands of 
Co** ions move —5 eV below Fermi levels, showed in the lower 
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Figure 2. Simulated X-ray diffraction for Li, ;CoSiO, and LiCoSiO,: Pmn2,-DP and Pbn2. 


panel of Figure 3. Because of the strong hybridization, a 
remarkable change in the Co(3d)-O(2p) orbital overlap (so- 
called rebonding shift) enhances the covalent bonding between 
the Co ions and the surrounding oxygen ligands. No obvious 
sharp occupied spin-down e-bands near the gap presents in the 
DOS of Co**, thus, it is a charge-transfer-type insulating. In 
LiCoSiO,, all the formal triple valence Co** ions are in the 
configuration e*t,* and the insulating gap becomes a complete 
charge transfer-type. The swapping of the interaction between 
3d(Co)- 2p(O) orbitals changes the nature of the insulating 
gap from Mott—Hubbard-type to charge-transfer-type, a similar 
phenomenon has been observed in other systems.“ This 
swapping is an evident proof of the band structure change by 
lithium intercalation/extraction, which is different from simple 
energy level shift in rigid-band model or redox couple pinning." 

The swapping of near-gap states correlates with the 
contraction of the oxidized CoO, units, establishing the 
cationic electrostatic polarization.'* The correlation can be 
elucidated through a process of tracking how the electronic 
structure changes along the structural relaxation after removing 
the lithium ions. Figure 4a shows that Pmn2,-DP DOS is a 
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Figure 4. Partial density of states (DOS) of delithiated: (a) 
Li, ;CoSiO,; (b) LiCoSiO,. The upper is for nonrelaxed structures, 
with the lower for well-relaxed structures and only shown for Pmn2,- 
DP. Fermi levels are at zero. 


metal-like before lattice relaxation, a clear gap-opening occurs 
after relaxing the structure of Li,;CoSiO,, along with the 
building-up of spin ordering by the oxidation of CoO, arrays. 
Figure 4b also shows the gap increasing after relaxing 
LiCoSiO,. Similar gap increasing presents in Pbn2,, too. For 
Li, ;CoSiO,, the gap of Pmn2,-DP opens from 0 to 1.3 eV, 
larger than the gap increase of Pbn2, from 0.6 to 1.5 eV; for 
LiCoSiO,, the gap increase of Pmn2,-DP from 0.6 to 1.5 eV is 
smaller than that of Pbn2, from 0.2 to 1.6 eV. 

The coexistence of both lattice distortion and band gap 
opening/increases demonstrates the occurring of Peierls 
distortions in both delithiated Li,;CoSiO, and LiCoSiO,, 
which are intrinsic Mott insulators. The resulting electronic 
localization with gap increases is the physical origin of difficult 
electronic transport during charging and discharging, which 
may account for the significant loss in discharging capacity of 
Co silicates. 


4. CONCLUSION 


In this work, we combine first-principles calculations with XRD 
simulations to systematically study the crystal structures, 
bonding characteristics and local electronic structures of 
Li,CoSiO, (x = 2.0, 1.5, 1.0) polymorphs with space groups 
Pmn2,-DP and Pbn2,. Average voltages by lithium extraction 
from Li,CoSiO,, through Li,;CoSiO, to LiCoSiO, are 
calculated, agreeing well with available experiment. The 
tetrahedron bonding network forms a special pattern in the 
lattices of Li,;CoSiO,, indicating the formation of spin 
ordering of CoO, arrays, and may be observable in XRD 
spectra according to simulation. Ordering of CoO, arrays can 
be used to characterize the redox process during lithium ions 
extraction/intercalation. Both delithiated Li,;CoSiO, and 
LiCoSiO, are intrinsic Mott insulators. The swapping of the 
interaction between 3d(Co)- 2p(O) orbitals changes electronic 
band structures from Mott—Hubbard-type to charge-transfer- 
type, which correlates with the contraction of the oxidized 
CoO, units, demonstrating the occurring of Peierls distortions 
that may be responsible for the significant loss in discharging 


capacity. 
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